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1. Introduction 

Dimension provides an indication of the size and complexity of a set and various 
kinds, such as box counting, packing and Hausdorff dimensions, play an important 
role in the study of fractals |50j . For example, the Hausdorff dimension of the Can- 
tor middle third set is log 2/ log 3 (proved by Hausdorff in his seminal paper |60j), 
that of the Koch snowflake curve is log 4/ log 3 and it has recently been shown that 
the boundary of the Mandelbrot set, a very complicated set of Lebesgue measure 
or null set in the complex plane with topological dimension 1, has Hausdorff di- 
mension 2 |122j . On the other hand, Diophantine approximation is a quantitative 
analysis of rational approximation and so, at least at first sight, is less geometrical. 
The purpose of this article is to show that Hausdorff dimension plays an important 
part in this theory too. 

In order to keep the article accessible, the emphasis is on approximation of real 
numbers by rationals and the less well known topic of approximation of complex 
numbers by ratios of Gaussian integers. The more general theory, which recently has 
seen some spectacular advances, will be referred to and some applications sketched. 
The article is organised as follows. We begin with a brief treatment of Hausdorff 
measure and Hausdorff dimension. We then explain some of the principal results in 
Diophantine approximation and the Hausdorff dimension of related sets, originating 
in the pioneering work of Vojtech Jarmk |98l . We conclude with some applications 
of these results to the metrical structure of exceptional sets associated with some 
famous problems. It is not intended that all the recent developments be covered 
but they can be found in the references cited. 



2000 Mathematics Subject Classification. Primary 11J83; Secondary 28A78, 37F50, 28A80. 

Key words and phrases. Diophantine approximation, exceptional sets, Hausdorff dimension, 
fractals, small divisors, dynamical systems, complex dynamics. 

Research funded by EPSRC grant no. GR/N02832/01 with additional support from INTAS 
grant no. 001-429. 



1 



©2003 American Mathematical Society 



2 



M. MAURICE DODSON AND SIMON KRISTENSEN 



2. Hausdorff measure and dimension 

Felix Hausdorff introduced the notion of Hausdorff dimension in a remarkable 
and influential paper .60, that extended Caratheodory's approach to Lebesgue mea- 
sure [24. in a simple but far-reaching way. (Kahane's Foreword to the book Frac- 
tals |29| includes a short and moving biography of Hausdorff.) Dimension had been 
taken to be a non-negative integer but by a simple observation, which Hausdorff 
described modestly as 'a small contribution', he modified Caratheodory's definition 
of measure to obtain a measure associated with a dimension that could be any non- 
negative real number. We shall assume a knowledge of Lebesgue measure and as 
usual, we shall often say almost no to indicate a null set - thus almost no numbers 
are rational - and we shall say almost all to indicate a set whose complement is 
null, so that almost all numbers are irrational. 

For familiar sets such as the interval, circle and the plane, the Hausdorff dimen- 
sion (defined below) coincides with the usual notion of dimension and is respectively 
1, 1 and 2. However, a significant difference is that any set in Euclidean space has 
a Hausdorff dimension (a non-measurable set in M" has full Hausdorff dimension 
n). In particular, null sets, such as Cantor's middle third set or the set of badly ap- 
proximate numbers (see t !3.2|) , have a Hausdorff dimension and this gives a way of 
discriminating between them. It is also natural to study the Hausdorff dimension of 
exceptional sets which are sets associated with the invalidity of some result, making 
it desirable that they be null. A brief and more or less self-contained account of 
Hausdorff measure and dimension is now given (more detailed expositions can be 
found in fl6ll49l f50l f52l f89l [T7)6| ) . 

2.1. Hausdorff measure. Caratheodory's approach to the measure of a set 
E in W 1 was based on 'approximating' E by countable covers consisting of small 
'simple' sets U in M". Hausdorff 's idea was to introduce for a given cover, 'rf say, 
of E the sum (sometimes termed the s-length of the cover c to) 

:= ^(diam[/) s , 

where diamC/ = sup{|x — y| 2 : x.y e U} is the diameter of U (|x — y| 2 is the usual 
Euclidean distance between x and y) and where s is a non-negative real number that 
is not necessarily an integer. Hausdorff considered a monotonic function I which 
allows more discrimination but for simplicity we shall stick to the more familiar 
widely used special case Z(diami7) = (diam[/) s , associated with what is now usually 
called the Hausdorff dimension but also sometimes called the Hausdorff-Besicovitch 
dimension |106j . The possibly infinite number t s (^?) gives an indication of the 
number of subsets U in ^ needed to cover E. In order to effect the approximation, 
the diameter of the sets U in the cover is restricted to be at most 6 > 0. 
Let 

H S S (E) := inf V (diam[/) s = in££ s (%), 

where the infimum is taken over all covers % of E by sets U with diam U ^ 5; such 
covers are called <5-covers. For a point x, 7i|({x}) = 1 when s = and vanishes 
when s > 0. As S decreases, 7i| can only increase as there are fewer ?7's available, 
i.e., if < S < 8', then 

H S S ,(E)^H S S (E). 
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The set function 7i s 5 is an outer measure on W 1 but the limit Tt s (which can be 
infinite) as 6 — > 0, given by 

(2.1) H S (E) = UmH s s (E) = sup H S S (E) e [0,oo], 

is better behaved. From its construction by covers, 7i s (E) ^ 7i s (F) for any E C F 
and indeed is subadditive and a regular outer measure. The restriction to the a- 
field of Tt s measurable sets (which includes open and closed sets, limsup and liminf 
sets and Gg and F a sets) is usually called the Hausdorff s- dimensional measure. 
Hausdorff 1-dimensional measure coincides with 1-dimensional Lebesgue measure 
and in higher dimensions, Hausdorff n-dimensional measure is comparable to n- 
dimensional Lebesgue measure, i.e., 

7i n (E) x \E\, 

where \E\ is the Lebesgue measure of E and where for a, b > 0, a x b means there 
exist constants c, d > such that a ^ cb ^ c'a or a = 0(b), b = 0(a) in Landau's 
O-notation. Thus a set of positive n-dimensional Lebesgue measure has positive 
Hausdorff n-measure. 

Because it is defined in terms of the diameter of the covering sets, Hausdorff 
s-measure is unchanged by restriction to closed, convex or open sets. It is also 
unchanged by isometries and so in particular by translations and rotations. It is, 
however, affected by scaling in the natural way (as are fractals): for any r > 0, 

H s (rE) =r s H s (E). 

2.2. Hausdorff dimension. Zero-dimensional Hausdorff measure Ti (E) is 
simply counting measure; thus the Hausdorff s-measure of a set of k points is k 
when s = and for s > 0. This pattern is typical. When the set E is infinite, 
7i s (E) is either or oo, except for possibly one value of s. To see this, the definition 
of W| (E) implies that there is a (5-cover % of E such that 

J2 (diamC) 8 ^ H%{E) + 1 < H S (E) + 1 < oo. 

Suppose that H s ° (E) is finite and s = sq + e, e > 0. Then for each member C of 
the cover (diamC) So+e ^ <5 e (diam C) s ° , so that the sum 

(diamC) 50 ^ < <5 e ^ (diamC) S0 . 

Hence 

H S S Q+£ {E) s$ J2 (diamC) S0+£ < 5 £ (diamC) 50 sC 5 e {H s «(E) + 1) 
and so 

U S (E) = H S ° +E (E) = \iiaH s 5 ° +£ (E) ^ lim S £ (H S °(E) + 1) = 0. 

On the other hand suppose H S °(E) > 0. If for any e > 0, H S °~ £ (E) were finite, 
then by the above H s o (E) = 0, a contradiction, whence H S ~ £ (E) = oo. 

Thus for each infinite set E in n-dimensional Euclidean space, there exists a 
unique non-negative exponent sq such that 



H S (E) = 




sC s < s , 

S < S < OO, 
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Figure 1. The graph of the Hausdorff measure of E. The Haus- 
dorff dimension dimn-E = sq, the point of discontinuity. 

as shown in Figure (reproduced with the permission of the Cambridge University 
Press from [16 ). The critical exponent 



where the Hausdorff s-measure crashes is called the Hausdorff dimension of the set 
E and is denoted by diniH-E- Thus the Hausdorff dimension of a finite set is 0, as 
it is for a countable set. It is clear that 

(2.3) if H S {E) = then dim H £ s; and if H S {E) > then dim H £ > s. 

The Hausdorff dimension tells us nothing about the Hausdorff s-measure at 
the critical exponent so = dim^E, only that this is the appropriate exponent to 
investigate the measure. The sudden change in Hausdorff s-measure at so = dimn^ 
can be compared to the focal length of a microscope. If the lens is too close, the 
image fills the eyepiece and cannot be resolved; if the lens is too far away, the image 
is invisible. At the focal length, the image is in focus and can be resolved. 

The main properties of Hausdorff dimension for sets in R n are 

(i) UECF then dim H £ < dim H F. 

(ii) divan E ^ n - 

(hi) If \E\ > 0, then dim H -E = n. 

(iv) The dimension of a point is 0. 

(v) If dirnn-E < n, then \E\ — (however dimni? = n does not imply \E\ > 0). 

(vi) dimn(-Ei x E-i) ^ dimH-Ei + dira^E-z 

(vii) dim H U™! Ej = sup{dim H £' :) : j e N}. 

It can be shown that the Hausdorff dimension of any countable set is and that 
of any open set in R™ is n ^Q, p. 29]. The nature of the construction of Hausdorff 
measure ensures that the Hausdorff dimension of a set is unchanged by an invertible 



(2.2) 



s = inf{s G [0, oo) : H S (E) = 0} 
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transformation which is bi-Lipschitz. This implies that for any set S C R \ {0}, 
dimuS 1-1 = dinin^, where S^ 1 = {s -1 : s e S}. Thus on the whole, Hausdorff 
dimension behaves as a dimension should, although the natural formula 

dim H (-Ei x E%) — dim H -Ei + dim H i?2 

does not always hold (49 §5.3] (it does hold for certain sets, e.g., cylinders, such 
as E x I, where I is an interval: dimn(£' x /) = dim^E + diniH-f = dimn-E + 1 by 
(iii), see [16]). 

The general character of i5-covers in the definition of Hausdorff outer measure 
can be difficult to work with and for many applications in higher dimensions, it is 
convenient to restrict the elements in the 5-covers of a set to simpler sets such as 
balls or cubes. For example, covers consisting of hypercubes 

H = {xeR n : |x-a|oo < 6}, 

where |x|oo := max{|x :) | : 1 < j ^ n} is the height of x G R", centred at a e E n and 
with sides of length 26 are used extensively. While outer measures corresponding 
to these more convenient restricted covers are not the same as Hausdorff measure, 
they are comparable and so have the same critical exponent 89, Chapter 5]. Thus 
there is no loss as far as dimension is concerned if the sets U are chosen to be balls 
or hypercubes. Of course, the two measures are identical for sets with Hausdorff 
s-measure which is either or 00. Such sets are said to obey a 'O-oo' law, this being 
the appropriate analogue of the more familiar '0-1' law in probability 76, p 339]. 
Sets which do not satisfy a O-oo law, i.e., sets which satisfy 

(2.4) <H dinmE {E) < 00, 

are called s-sets; these occur surprisingly often |501 p. 29] and enjoy some nice 
properties (see Chapters 2-4 of [49j). One example is the Cantor set which has 
Hausdorff s-measure 1 when s = log2/log3 49, p. 14]. However it seems that 
s-sets are of less interest in Diophantine approximation where the sets that arise 
naturally, such as the set of badly approximable numbers or the set of numbers 
approximable to a given order (see next section), obey a O-oo law. The first steps 
in this direction were taken by Jarmk, who proved that the Hausdorff s-measure of 
set of numbers rationally approximable to order v (see M.2I1 was or 00 |67L I68j . 
This result turns on an idea related to density of Hausdorff measure. 

Lemma 2.1. Let E be a null set in R. Suppose that for any interval (a, b) and 
« 6 [0,1], 

(2.5) H s (ED(a,b)) ^ K {b - a)H s (E) . 
Then H S (E) = or 00. 

Proof. Suppose the contrary, i.e., suppose < H S (E) < 00. Since E is null, 
given e > 0, there exists a cover of E by open intervals (cij, bj) such that 

Y^(bj -aj)< e. 

3 

By H2.5H , there exists a constant K > such that 

< H s (E) = H s (U j (aj , bj ) n E) < KH S (E) ( b j - «j ) < K ^ s (E) <H S (E) 

3 

for e < l/K, a contradiction. □ 
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The proof for a general outer measure is essentially the same. The sets we 
encounter in Diophantine approximation are generally not s-sets and some satisfy 
this 'quasi-independence' property. 

For other definitions of dimension, such as box-counting and packing dimension, 
and their relationship with Hausdorff dimension, see ,50, Chapter 3]. 

2.3. The determination of Hausdorff dimension. Unless some general 
result is available, the Hausdorff dimension dinin^ of a null set E is usually deter- 
mined in two steps, with the correct upward inequality Aiva^E ^ so and downward 
inequality dimn-E ^ so being established separately. 

2.3.1. The upper bound. In view of 112. 3|) . an upper bound can be obtained by finding 
a value of s for which H S (E) vanishes. To find such a value, it suffices to exhibit 
a cover {H} of E (E C Off) by hypercubes H of arbitrarily small sidelength and 
s-length. This can often be done by adapting the estimate involved in showing that 
Lebesgue measure is 0. When E is a limsup set, i.e., 



oo oc 



E = lim sup E N = P| [J E k = {x £ K™ : x £ E k for infinitely many k £ N} 

N ^°° N=lk=N 

for a sequence of sets E n , a simple Hausdorff measure counterpart of the convergence 
case of the Borel-Cantelli lemma often gives the correct upper bound for dimnE. 
This is useful in Diophantine approximation. 

Lemma 2.2. Let 

E = {x £ R n : x £ Ek for infinitely many k £ N}. 

If for some s > 0, 

oo 

(2.6) ^diam(E*) s < oo, 

fc=i 

then H S (E) = and dim H £ < s. 

Proof. From the definition, for each N = 1,2,..., 



£cy E k , 

k=N 

so that the family = {E k : k ^ N} is a cover for E. By iETfij) . 

OO 

lim diam(E k ) s = 0- 

k=N 

Hence lim^oo diam(-E'fc) = and therefore given 5 > 0, ff( N ^ is a <5-cover of E for 
N sufficiently large. But 

oo 

Hg(E) = inf ^ (diamcT) 5 ^ l s { c g {N) ) = ^ diam(£; fc ) s -> 
as N -> oo. Thus H S 5 (E) = and by = 0, whence dim H £ s. □ 



HAUSDORFF DIMENSION AND DIOPHANTINE APPROXIMATION 



7 



2.3.2. The lower bound. The lower bound is often harder (though by no means 
always, see 39|) . It requires showing that given any s < sq and any cover {C} 
of E with the diameters of the covering elements arbitrarily small, the s-length 
^ c (diamC) s ^ 5 for some positive 6. This can be very difficult and has led 
to the development of a variety of methods. In Diophantine approximation, the 
regular systems introduced by Baker and W. M. Schmidt jjjjj and the more general 
ubiquitous systems [45\ depend on a good supply of approximating elements (e.g., 
the rationals) . These and related techniques have proved effective in obtaining lower 
bounds for the Hausdorff dimension of sets of number theoretic interest (see the 
survey article |13j for more details) . A more fundamental approach is the so-called 
mass distribution principle. 

Lemma 2.3. Let [i be a measure supported on a bounded Borel set E in W 1 . 
Suppose that for some s ^ 0, there are strictly positive constants c and S such that 
H{B) < c (diamB) s for any ball B in W l with diamB 6. Then H S {E) ^ fi(E)/c. 

PROOF. The proof is short. Let {B k } be a (5-cover of E by balls B k . Then 

fi(E) < (i l\jB k J <^ M (^fc) ^c^(diamB) s . 
V k ) k k 

Taking infima over all such covers, we see that H$(E) ^ fi(E)/c, whence on letting 

H S (E) ^n(E)/c>Q. 

□ 

This simple lemma is surprisingly useful and gives the easy part of Frostman's 
lemma [56] which is now stated in full. The Vinogradov notation a -C b for a, b > 
means that a = 0(b). 

Lemma 2.4. Let E be a Borel subset ofW 1 . Then 

H S (E) > 

if and only if there exists a measure [i on K™ supported on E with fJ,(E) finite such 
that /i(B) <§C (diamS)" for all sufficiently small balls B. 

Thus if E supports a probability measure \i (fi(E) — 1) with [i(B) (diam£?) s 
for all sufficiently small balls B, then dimn-E ^ s. The converse is more difficult 
but can be proved using net measures (see |25L I50L [89 i. 

3. Diophantine approximation 

At its simplest level, Diophantine approximation is concerned with approximat- 
ing real numbers by rationals. Hardy and Wright's classic Introduction to the theory 
of numbers 57j contains an excellent account while the more advanced (261 H19j 
are devoted wholly to Diophantine approximation. The theory extends to approx- 
imating vectors in R™ (simultaneous Diophantine approximation) and to matrices 
(systems of linear forms). For simplicity, we will stick mainly to one particular 
direction in the one dimensional real and complex cases and treat the extensions 
to higher dimensions and other settings fairly briefly. Since the rationals Q are a 
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dense subset of the real numbers R, given any real number a and any positive e, 
there exists a rational p/q such that 



(3.1) 



< e. 



The numerator p is often of no interest and the size of the expression 
(3.2) \\qa\\=zmn{\qa-p\:peZ}, 

the distance of qa from the integers Z, is considered. Although convenient, it will 
not be used much here in order to keep the notational burden to a minimum. 

In simultaneous Diophantine approximation, one considers the system of n 
inequalities 

Pk 



oik 



Q 



< e, k =, 1 



This system can be expressed more concisely as a single vector inequality with 
a = (ai, . . . , a n ) € R™, p € Z n , q € N, by considering the expression 

a-E 

q 

where for x £ R™, = max{|a;i|, . . . , |x n |} is the height of x, or, on multiplying 
by q, the expression 

\qa - p|oo = \\qa\\ = max{||ga,-|| : j = 1, . . . , n}. 

The last inequality has a dual or linear form version: given a £ R™ and e > 0, one 
considers the inequality 

|q-a-p| <£, 

where q € Z" and p e Z. The last two inequalities can be combined into a single 
general one. The system of n real linear forms 



J 



1. 



in to real variables £i, . . . ,£ n , can be written more concisely as £A, where A — (a^) 
and the system of n inequalities in to variables 



\qiaii 



+ q m a m \ -Pl\ < £ 



\qiai„ H h q m a mn — p n \ < e i 

can be written \<\A — p|oo < £• Further details are in [26l l57l ITT9] . The theory 
extends naturally to the fields of p-adic numbers [27|, 18 7| and formal power se- 
ries |82U84| . Less obviously, it also extends to discrete groups acting on hyperbolic 
space. This is relevant to Diophantine approximation over the Gaussian integers or 
rationals considered below in *14.3I so an outline is now given. 

The hyperbolic space setting sprang from the observations that the the real axis 
is the set of limit points of the rationals and that the rationals can be characterised 
as the parabolic vertices of the modular group SL(2,Z), i.e., as the orbit of the 
point at infinity under the linear fractional or Mobius transformations 

(3.3) — a,b,c,d£Z, ad—bc=l, 
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of the extended upper half plane H 2 ~ {z ~ x + iy e C: y ^ 0} U {oo}. For each 
element 



in SL(2, Z), the point oo in the extended real line evidently goes to a/c under the 
group action and is also a fixed point of the map z i— > z + 1. The maps g form 
the modular group, SL(2, Z), a discrete subgroup of SL(2,R), which is essentially 
(modulo the centre) the group of orientation preserving Mobius transformations 
of the upper half plane H 2 to itself. When the upper half plane is endowed with 
the hyperbolic metric derived from dp — \dz \ /y, it is a model for two dimensional 
hyperbolic space (H 2 ,p). The Mobius group M(H 2 ) is the group of isometries of 
(H 2 ,/?). Because the group SL(2,Z) is discrete, points in the orbit can accumulate 
only on the boundary R U {oo} of H 2 and because the group elements are isometries 
with respect to the hyperbolic metric, the limit set of any orbit is the extended real 
line R U {oo}. A discrete subgroup of M(H 2 ) is called a Fuchsian group. Further 
details of this rich and beautiful theory are in O |4l I11L 1951 [99j and there is a 
short survey in Chapter 7 of [l6j. 

These observations allow the classical theory, including the metrical theory, of 
Diophantine approximation to be translated into Fuchsian groups acting on the 
hyperbolic plane and to the much more general setting of Kleinian groups acting 
on (n + l)-dimensional hyperbolic space (M n+1 ,p), n ^ 2 (Kleinian groups are the 
discrete subgroups of the Mobius group M (H™ +1 ) of isometries of H n+1 ; further 
details are in the references given). The Picard group SL(2,Z[i]) consists of 2 x 2 
matrices over Z[i] with determinant 1 and has an action on C given by 

(3.4) z i ► — a ' ^' c ' ^ e ^M' o-d — bc— 1. 

The limit set of the Picard group is the extended complex plane (or Riemann sphere) 
C U {oo} and the orbit of the point at oo under the group is the set of ratios of 
Gaussian integers [9lj . Thus the Picard group plays a role precisely parallel to that 
played by the modular group, expressed by in approximating real numbers 

by ratios of integers. 

In the literature cited, hyperbolic space is usually taken in the equivalent 
Poincare form of the open unit ball B ri+1 = {x G W 1+1 : |x| 2 < 1}, where the ball is 
now endowed with the equivalent hyperbolic metric p given by dp = \dx\ 2 /(I— l*^)- 
To ease comparison, we shall adopt this viewpoint, even though the upper half plane 
model, as used by Sullivan in |130j . is more natural for Diophantine approxima- 
tion. We choose to consider the (n+ l)-dimensional hyperbolic space, as the results 
in Diophantine approximation are results about the boundary of H 11+1 , which is 
rt-dimensional. 

The analogue p of the point at infinity for Kleinian groups is not quite straight- 
forward. First of all the nature of the elements g of the Kleinian group G implies 
that each g has at most two fixed points on the boundary of the ball. The special 
point p is called a parabolic fixed point if it is the unique fixed point on the bound- 
ary of some element in G; otherwise they are called hyperbolic fixed points. The 
orbit of a special point p under the action of a Kleinian group G corresponds to the 
rationals Q. The limit set A(G) of the orbit under G of a point in H™ +1 lies in the 
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boundary S n , Given a £ A(G) C §™, one considers the quantity 

\<*-g(p)\2, 

where | ■ |a is the usual Euclidean metric in K n+ , Analogues of the principal 
theorems in Diophantine approximation have been obtained with relatively minor 
technical restrictions and will be discussed below; a brief survey is in Chapter 7 
of |16| . There is a striking dynamical interpretation of the approximation in terms 
of flows on the associated quotient space H n+1 /G; more details are in |16LI53ll9ll 
11301 IT35] . We now return to the one-dimensional theory. 

3.1. Dirichlet's theorem. It is not difficult to make 13. If) more precise: given 
any real number a and any positive integer q, there exists an integer p such that 
\qa — pj < 1, and indeed such that 



P 

a 

q 



1 



There are denominators q for which more can be said by using Dirichlet's celebrated 
'box argument' (see 1261 157j l. 

Theorem 3.1. For each real number a and any positive integer N > 1, there 
exists a rational p/q with denominator satisfying 1 ^ q ^ N, such that 



P 



1 1 

< ~qN ^ <f 



Proof. Let [a] be the integer part of a and {a} its fractional part, so that 
a = {a] + {a}. Divide the interval [0, 1) into N subintervals [k/N, (k + l)/N), where 
k = 0, 1, . . . , N-l, of length l/N. The N+l numbers {ra}, r = Q,l,...,N, fall into 
the interval [0,1) and so two, {ra}, {r'a} say, must fall into the same subinterval, 
[k/N, (k + 1)/N) say. Suppose that r > r'. Then 

|{ra} — {r'a}| = |ra — [ra] — r'a + [r'a] \ = \qa — p\ < — , 

where q — r — r',p — [ra] — [r'a] G Z and 1 ^ q ^ N. Dividing by q gives the 
quantitative inequality 



(3.5) 



< 



qN 



and since 1 ^ q ^ N, the final inequality is immediate. 



□ 



A nice sharpening is in [58, p. 1]. When p, q are restricted to having highest 
common factor 1, the inequality 



(3-6) 



< 



1 



has only finitely many solutions if and only if a is a rational. Thus almost all 
real numbers satisfy 113. 6JI for infinitely many rationals p/q. Without this restric- 
tion, 13. (ill holds infinitely often for all a £ R. 

Dirichlet's theorem is one of the fundamental results in the theory of Diophan- 
tine approximation. It can be viewed as a result about covers and plays a central 
part in the Jarmk-Besicovitch theorem, discussed below. The theorem generalises 
to the simultaneous Diophantine approximation of n real numbers a±, . . . ,a n |57l 
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Theorem 200] and asserts that given N G N, there exists geff, with q < N and 
p = (pi, . .. ,p n ) G Z n such that 



P 

a 



< 



1 



In particular, for simultaneous Diophantine approximation in the plane, given any 
a = (ai, as) G M 2 and A^ G N, there exists q G N with q ^ N and Pi,p2 € Z such 
that 



(3.7) 



Pi 

ai 

q 



P2 

Oi% 

q 



< 



1 



qN 1 / 2 



-3/2 



There is a so-called dual version: given a G R™ and N G N, there exists q G Z ra 
and a p G Z such that 

|q-a-p| < iV"™. 

These can be combined into a result for systems of n real linear forms ; 221 Chapter 1, 
Theorem VI]: 

\qA - < N~ m ' n . 

In the setting of a Kleinian group acting on hyperbolic space, the analogue of 
the denominator in Theorem 13 . II corresponding to g{p) is defined to be 



(3.8) 



A 9 := Idct^lo)!- 1 = icoshp(0, 5 (0))xe^°)) 



in the ball model, and so A g — > 00 as |g(0)| — > 1, i.e., as the orbit of the origin 
moves towards the boundary. Here Dg\o denotes the Jacobian of g evaluated at the 
origin. For finitely generated Fuchsian groups of the first kind taken to be acting 
on the closed unit disc A, the elements g are of the form 



9 



- b I, a,beAM-\br = l 



and A s = 2(|a| 2 + |6| 2 ). Hedlund's lemma, which is a classical result in the theory 
of discontinuous groups, is a partial analogue of Dirichlet's theorem. Let ( £ S 1 , 
the unit circle. For any £ G S 1 which is not a parabolic point, there exist infinitely 
many g G G such that 

1^-9(01 <f , 

for some C > 0. A complete analogue of Dirichlet's theorem, including the quan- 
titative inequality l|3,5Jl . was obtained by Patterson [99] for Fuchsian groups and 
later he and others extended it to Kleinian groups [lflfll 11281 1129L fl35| . However, 
the statements in the Kleinian group setting differ for parabolic and hyperbolic 
fixed points and so for simplicity the result will be stated when G has a unique 
parabolic point p. Let N ^ 2. Then for any £ G A(G), there exists a g G G with 
X g < N such that 



l£-<?(f>)| 2 < 



C 



where C is a constant depending only on G. 
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3.2. Types of approximation. The equation H3.6JI given by Dirichlet's the- 
orem is essentially best possible as Hurwitz [64\ showed that each a £ K satisfied 
the inequality 



P 

a 

q 



< 



TV 



for infinitely many positive integers q and that this was best possible in the sense 
that the constant 1 / V5 cannot be reduced for numbers a equivalent to the golden 

ratio (V5 + 1)/2 f26i[55inrrirrT9| . 

3.2.1. Badly approximable numbers. A number j3 € M is called badly approximable 
or of constant type if there exists a K — K(j3) such that 



(3.9) 



q 



K 



forallp/g € Q. Using H3.2I) . H3.9I) ca,n he written as a\\a3\\ ^ K . In view of Hurwitz 's 
theorem, the constant K < l/v5- Quadratic irrationals, such as v2 and the golden 
ratio (v5+ l)/2, are badly approximable. This is proved in |57l §11.4] using the 
fact that the partial quotients in the continued fraction expansion for a quadratic 
irrational are periodic. However, the proof relies only on the boundedness of the 
partial quotients, which therefore characterises the badly approximable numbers. 
The set of badly approximable numbers will be denoted by <B. The notion extends 
naturally to higher dimensions and to the more general settings mentioned above. 

Badly approximable numbers are important in applications, particularly in sta- 
bility questions for certain dynamical systems |121| . For example the 'noble' num- 
bers, which are equivalent to the golden ratio, have been conjectured to be the 
most robust in the breaking up of invariant tori [88]. One very practical applica- 
tion involved the design of rocket casings. These were made using ruled surfaces 
and vibrations from the motors were propagated along the generators. To reduce 
the effects of resonance and delay the onset of catastrophic vibration, the ratio of 
the circumference to the length of the casing was chosen to be a quadratic irra- 
tional (V. I. Arno'ld, personal communication). The desirable properties of badly 
approximable numbers (and in particular of the golden ratio) appear to be related 
to their occurrence in nature. It has recently been discovered that the ratio between 
two step heights on the surface of certain quasi-crystals is given by the golden ratio 
(see |3fl| for statements of this result and additional examples) . 

The notion of badly approximable numbers carries over to higher dimensions, 
including systems of linear forms |118j , p-adics [jQ and fields of formal power series 
82j as well as to Kleinian groups acting on hyperbolic space 20 [ I53L I99L llOl] . In 
the hyperbolic space setting, a point (5 in A(G) is said to be badly approximable 
with respect to p if there exists a positive constant K — K((3) such that 

\P-g(p)\ 2 >K/\ g 

for all g € G. 

3.2.2. Diophantine type. The concept of a badly approximable number has exten- 
sions to restricted classes of real numbers and points in R™ that are useful in 
connection with stability and other questions (see © and fortunately enjoys full 
measure. Let K > 0, v > 1. The real number a is said to be of Diophantine type 
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(K, v) [gj (the definition has been altered slightly for consistency) if 

K 



for all rationals p/q; the set of numbers of Diophantine type (K, v) is denoted by 
T>(K,v). The union 

V v = \J V(K,v) 

K>0 

consists of numbers of Diophantine type v (i.e., of type (K, v) for some K > 0) and 
the union 

V=\JV V = |J V(K,v) 

v>l K>0,v>l 

is the set of numbers of Diophantine type v for some v > 1. Note that X>i = 05, 
the set of badly approximable numbers. We will see in iQ]that 2} is null but that 
when v > 1, T> v has full measure, which is pleasing since points of Diophantine 
type have desirable approximation properties for certain applications. Again, the 
notion of Diophantine type extends naturally to higher dimensions |46j . We will be 
particularly interested in the case of a single linear form and accordingly we extend 
the definition of Diophantine type for a real number to a point in R n . A point 
(3 € R n is of (dual) Diophantine type (K, v) if for all p € Z and non-zero q€ Z™, 

K 

q-/3-p > — — . 

|q|So 

3.2.3. Well approximable numbers. In applications, we will be interested in points 
a which are not of Diophantine type (K, v) for any K > 0, i.e., in one dimension 
with the set 
(3.10) 



a e R: for any K > 0, 



< 



K 



tor some - € 



= R\V V , 



the complement of V v . These numbers are closely related to numbers which are 

rationally approximable to order v+1 |57l §11.4], i.e., to the set 

(3.11) 



Rv 



a € R: for some K > 0, 



K r ■ n ■ , P 

< for infinitely many — 6 



(the exponent u + 1 in the denominator is a normalisation to keep the notation 
the same as elsewhere.) By Dirichlet's theorem, all real numbers are rationally 
approximable to order 2 and quadratic irrationals are rationally approximable to 
order exactly 2. 

The constant K in the definitions of R v and E v is of less significance than the 
exponent, which motivates the next definition. A number a which satisfies 

1 



(3.12) 



< 



q L 



for infinitely many p/q £ Q will be called v- approximable; if v > 1, a is called very 
well approximable. Thus Liouville numbers, which satisfy l|3.12H for any v, are very 
well approximable. We let W v denote the limsup set of w-approximable numbers, 
i.e., 

p 



W v = 



ja € R: 


P 


a 




q 



< v+1 for infinitely many — G 
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The sets E v , R Vl W v are related and decrease as v increases, as follows from 
the inclusions 



(3.13) 



Rv+e c E v c W v c R v , 



where e > is arbitrary. Clearly W v C R v . Next consider the complementary 
inclusion K \ W v C V v and let a W v . Then for all but finitely many positive 
integers q, 

j = l,...,r. 



P 

a 

q 



nV + 1 ' 



Moreover, for the exceptional values of q, say 



7 (r) 



K 



lin | 



q(i) a _ pii) 



(q U) ) V ■ i = l,...,r} >0, 



since qa £ Z implies fcga € Z for each fc 6 Z. Thus |a — ^ K/q v+1 for all 
p/<7 G Q and so a € V v , i.e., a ^ 

To establish the final inclusion, let a £ R v +e, so that for some K > 0, 



P 

a 

q 



< 



K 



q 



v+l+e 



for infinitely many p/q £ Q. Given any K' > 0, choose as we may, a denominator 
<7o sufficiently large so that q^ E ^ K'/K. Then there exists a po & Z such that 

K' 



Po 

qo 



< 



9o +1 



and a £ E v . The Jarmk-Besicovitch theorem (23 gives us the Hausdorff dimen- 
sion of W v and allows us to deduce from l|3.13fl that all the sets have the same 
Hausdorff dimension. 

3.2.4. ^ -approximable numbers. Now we look at the set of points which enjoy a 
more general approximation by rationals. A function : N — > M + such that 

lim *(q) = 

q — >oo 

will be called an approximation function; without loss of generality we can take 
^(q) ^ l/(2g) and later we shall also assume that q^(q) is decreasing (by decreasing 
we mean non- increasing here and subsequently). A real number a is said to be 'I' - 
approximable if a satisfies the inequality 



(3.14) 



P 

a 

q 



<*(<?) 



for infinitely many p/q € Q (there should be no confusion with v-approximable 
numbers defined above) . Note that this should not be confused with the inequality 
\qa — p\ < ip(q) which is often considered, particularly in higher dimensions and 
that there are other definitions depending on the form of H3.14I) . see for example 
|17| . The set of ^-approximable numbers in E will be denoted W(^). 
Since for each k € Z, (p + kq)/q £ Q and 



k- 



p + kq 





p 




a 




q 
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it follows that W(V) n[k,k+l) = {W(V) [0,1)) + k, so that W(V) can be 
decomposed into a union over unit intervals: 

= |J (W(¥) n [fe, jfc + i)) = |J n [o, i)) + k. 

k£l fcGZ 

As with the other types of approximation, the definitions of very well approx- 
imable, w-approximable and ^-approximable numbers extend naturally to systems 
of linear forms, p-adic numbers, formal power series and to the hyperbolic setting. 
We will be interested in the case of a single linear form: we say with some abuse of 
notation that a vector a. £ R" is w-approximable if 

(3-15) \q-a-p\<j-±- 

holds for infinitely many q G Z n and p £ Z; a is very well approximable if v > n. 
In the hyperbolic space setting, a point a in I" is ^-approximable with respect 
to G and p if 

l«-$(j>)la < *( A s) 

for infinitely many g 6 G. These and the other analogues of the real one-dimensional 
case will be discussed further when they arise below. 

4. Khintchine's theorem and metrical Diophantine approximation 

Following earlier work of Borel |2JJ , Khintchine gave an almost complete answer 
to the solubility of 113.1 4|) in terms of the measure of W(^). In a series of papers in 
the 1920's on the Lebesgue measure of the sets and <8 (22 EE EH EH , he 

laid the foundations of metrical Diophantine approximation. This theory, which is 
closely related to probability, measure and ergodic theory, considers sets of solutions 
to Diophantine inequalities in terms of Lebesgue and other measures. As a result, 
0-1 laws are a feature of the Lebesgue part of the theory, as in Khintchine's theorem 
below (see also |581 §2.2]). In addition, because an exceptional set for which a result 
is invalid can be of measure zero, this can lead to theorems having a strikingly simple 
yet general character. 

Theorem 4.1 (Khintchine). 

Jo, if Er=i < °°. 



|W(*)n[o,i) 



1, if fc^'(fc) is decreasing and Ylk=i k^{k) 



Thus M /r ( v l') is null when the series ^ fe k^>{k) converges and is full when k^(k) 
decreases and the sum diverges. Subsequently Khintchine extended the result to 
simultaneous Diophantine approximation [75 (see also [26 ) and Groshev extended 
it to systems of linear forms |126| . In particular the measure of the set of points 
(ai, a%) £ [0, l) 2 such that 



max 



Pi 

Q'l 

q 



P2 

a 2 

q 



<*(<?) 



for infinitely many pi,P2 £ Z, q £ N is or 1 accordingly as ^ k 2 ^(k) 2 converges 
or as fcvt'(fc) decreases and the sum diverges (c/. Theorem 14.61 below) . 

The convergence case follows readily from the Borel-Cantelli lemma |if6] and 
the argument is closely related to that of Lemma l2~2l However, the case of diver- 
gence is much more difficult and relies on a crucial 'pairwise quasi-independence' 
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result combined with 'mean and variance' ideas or with results from ergodic theory. 
For further details of this remarkable improvement of the Borel-Cantelli lemma, 
see [3ll l32ll58lFf9l Il05l Il26l 1131] (the books [HI Il26j also include accounts of 
W. M. Schmidt's important quantitative extension of Khintchine's theorem 116J). 
In terms of probability, we recall that the events Ej, j = 1, . . . , oo, with correspond- 
ing probabilities P{Ej) are pairwise independent if for any j ^ k, the probability 
of the two events occurring is given by 



As a result, the divergence half of the Borel-Cantelli Lemma holds for pairwise 
independence. Total independence, where the probability of any finite sequence 
of events is given by the product of the individual probabilities, is not necessary. 
Recently it has been shown that the divergence case is also related to the lower 
bound for Hausdorff dimension (see . 

There has been dramatic progress in the metrical theory over the last decades, 
particularly in the theory of 'dependent variables' where the point a lies on a 
manifold, so that coordinates of the point are functionally related. Sprindzuk's 
solution |127| of Mahler's conjecture in transcendence theory in terms of the Dio- 
phantine approximation of points on the Veronese curve {(x, x 2 , . . . , x n ) : x S K} 
gave this topic an enormous impetus which has seen the recent proof of Sprindzuk's 
conjectures and further results |12L I17L FT8j. 

The p-adic analogue of Khintchine's theorem was obtained by Jarnlk [751 an d 
extended to systems of linear forms by Lutz |87| . For fields of formal power series, 
the analogue was obtained by de Mathan |37j and has recently been extended 
to systems of linear forms [84]. The complex analogue of Khintchine's theorem, 
in which the approximation is by ratios of Gaussian integers, is discussed in H4.3I 
below. 

Khintchine's theorem corresponds to our intuition since if the approximation 
function \I> is large, then there is a better chance of the inequality being satisfied. 
In particular, the set W v is null for v > 1, since the series J2k converges, and 
full for v ^ 1 (when W v = [0, 1) by Dirichlet's theorem). 

Corollary 4.2. 



The result for the other sets follows from (13. 1311 . Less obviously, the theorem 
shows the Lebesgue measure of the set of a € [0, 1) such that H3.14J1 has infinitely 
many solutions is 1 when ^(k) = l/(fc 2 logfc) and when #(k) = l/(fc 2 (log k) 1+£ ) 
for any positive e. 

As a consequence of Jarnik's theorem on simultaneous Diophantine approxima- 
tion to be discussed in ii4.2l none of the sets W v , R v , E v is an s-set (see 112.411 for the 
definition). However, using Lemma T2. II and invariance under rational translates of 
R v , Jarmk had shown earlier that the set R v of numbers rationally approximable 
to order v is not an s-set and obeys a '0-oo' law |67L I68| . Although superseded by 
the above theorem, the argument is very nice but the papers cited are not read- 
ily available, so the proof is repeated here. Let k be a positive integer and let 



P{E j f\E k ) = P{E j )P{E k ). 



|w„n [0,1)1 = \Run [o,i)| = |K n [o,i)| 




when v > 1, 
when v ^ 1. 
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a e R v n [0, 1/fc], so that for some K = K(a) > 0, 



< 



K 

-,v+l 



for infinitely many rationals p/q. For each j/k, 1 ^ j ' ^ k — 1, 









a 




a 


9 







l_p_l 

k q k 



< 



Kk v+1 



{kq)^ 1 ' 

Thus a + j/k G R V D [j/k, (j + l)/k] and it follows that for each j = 0, 1, . 



,fc-l. 



iL n 



i i±i 

fc' fe 



= i?„ n 



fc' 



Moreover since Hausdorff measure is translation invariant, 

k-l 

/ : ) IN 

= feW 8 R v n 



H s (i^n [0,1]) =53^ s (-R„n 

whence for any j, k with 1 ^ j ^ k — 1, 
H " ( Rv n I" ' 1 ' 



j i + i 



fe' 



= -t;H s n [0, 1])) . 

Now every open interval (a,b) can be represented as a union of a countable set of 
intervals [j/k, (j + l)/fc], i.e., 



(o,6) = U 



jfc' fc 



that b — a = ■ fe l/fe. Hence 



H'((a,b))=H' (U 



i i + i 
fe' fc 



and since 7i s (-) is an outer measure, 

n s {(a, b)nR v n [o, 1]) = n s I |J R v n 



v j,fc 



j j + i 

fe' fc 



n 



j_ j + l 

fc' fe 



< n s (R v n [o, l]) 2 r = ( 6 - n [o, 1]). 



Thus the hypotheses of Lemma l2.1l are satisfied and so H s (R v r\[0, 1]) = TL S {R V ) is 
or oo. The argument can be extended to show that the set of numbers approximable 
to order v by algebraic irrationals is not an s-set. The measure at the critical 
exponent was shown to be oo by Bugeaud [23] . 

Khintchine's theorem also implies that the set 2$ of badly approximable num- 
bers is null. For given any K > 0, the sum J2 q (K/q) diverges and so by Khintchine's 
theorem the set of real numbers a satisfying [a — p/q[ < K/q 2 for infinitely many 
p/q € Q is full. Thus the complementary set F(K) of the set of a such that 
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\a — p/q\ ^ K/q 2 for all but finitely many p/q is null and evidently increases as K 
decreases. Prom its definition, 



35 C |J F(K) = |J F(l/N), 



K>0 N=l 

a countable union of null sets, whence 03 is null. 

Since the set W v of very well approximable numbers is null, the inclusions i|3.13Jl 
imply that V v is full for v > 1, i.e., that almost all real numbers are of Diophantine 
type (K, v) for some positive K. Thus almost all numbers are neither well or badly 
approximable. It is a remarkable fact, proved by Jarmk [65j in 1928, the same 
year as Besicovitch's first paper on Hausdorff measure and dimension (on planar 1- 
sets) |18j . that although the set 03 is null, dimuOS = 1, i.e., its Hausdorff dimension 
is maximal in the sense that it coincides with that of the ambient space M. This 
we now discuss. 

4.1. Jarrrik's theorem for badly approximable numbers. Let S (0, 1) 

and let a n , n — 1,2,..., denote the partial quotients of the continued fraction 
expansion for 9. For each N &N, define 

M N := {6: a n < N}. 

Now a number is badly approximable if and only if it has bounded partial quotients 
a n [57\ §11.4], so that 23 = limN^oo Mjy. In a pioneering paper that was the first 
on Hausdorff dimension in Diophantine approximation, Jarmk [55^ showed that for 
each N ^ 8, 

4 1 
1 _ TFi — 7i ^ dim H MAr < 1 



7Vlog2 a " 8iVlog7V 

It is of course immediate that dimn23 ^ 1 since 03 C K. We now state Jarnlk's 
theorem for badly approximable numbers, which follows from the above. 

Theorem 4.3 (Jarmk). 

dim H 03 = 1. 

Using (a, (3) games, W. M. Schmidt |117L ITlS] proved much more, extending 
Jarnlk's theorem to higher dimensions, so that the set of simultaneously badly 
approximable points in the plane has Hausdorff dimension 2. In fact, 03 is a thick 
set. This is a 'local' property in the sense that for each open interval /, 03 has 'full' 
Hausdorff dimension, i.e., dimnOSn/ = 1. A very general inhomogeneous analogue 
has been proved using quite different ideas from dynamical systems |77| . 

The game in one dimension involves two players A and B, a non-empty set 
Set and two parameters a € (0, 1), given to the player A, and f3 G (0, 1), given 
to the player B. Player B begins by picking a closed interval B\. Then A chooses a 
closed subinterval A\ C B\ with \A\ \ — a\B\\. Then B picks an interval B 2 C A\ 
with \Bi\ = f3\A\\ — f)a\B\\ and then A chooses another subinterval A 2 C B 2 with 
\A 2 1 = a\B 2 \ = a 2 /3\Bi \ and so on. Clearly the intervals B\, A\,B 2 , A 2 , . . . , form a 
decreasing nested sequence so that their intersection is a point, to say, in B\. Player 
A is called the winner if DjAj = {to} C S, otherwise B wins. 

A fuller account of the game in C will be given in so we will simply say that 
Schmidt showed that when S = 03, A can force w to be badly approximable, even 
though 03 is null, and deduced that dimH03 ^ 1. Since 03 C K, dimnOS = 1. Note 
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that since Ti}(*&) = |58|, the Lebesgue measure of 25 and since |23| = (Khintchine's 
theorem), it follows that the Hausdorff measure of 58 vanishes at s = diniH^B and 

W s («8) = 

Thus 55 is not an s-set (this can also be proved using Jarnik's lemma in t !2.2l above). 

Jarnik's theorem has also been extended to p-adic fields QJ, to fields of formal 
power series |82| and to hyperbolic space. It follows from the hyperbolic space 
counterpart of Khintchine's theorem that for geometrically finite Kleinian groups 
G, the set 25(G, p) of the hyperbolic analogue of badly approximable points, has 
zero Patterson measure and the analogue of Jarnik's theorem holds |53L 1991 llOlj . 
i.e., 

dim H s B(G,p) = dim H A(G). 
In addition, the exponent of convergence of G, 

5(G) := mf{s > 0: ^ \~ s < oo} = dim H A(G) 

seG 

|20L I100L In view of Lemma l2~2l this is perhaps not so surprising. In a 

striking parallel with continued fractions, badly approximable points correspond to 
bounded orbits of flows on manifolds [20l [3H QOH [53l Il20] . 

4.2. Jarrrik Besicovitch theorem. When v > 1, the set T> v is complemen- 
tary to the set R v of points approximable to exponent v, which is related to the 
set W v of w-approximable numbers (these statements also hold for the higher di- 
mensional analogues) . The Hausdorff dimension of W v was determined by Jarnik 
in 1929 [66| and independently by Besicovitch in 1934 (19] . 

Theorem 4.4 ( Jarnlk-Besicovitch) . When v ^ 1, 

2 

dmin^ = — — , 
v + 1 

and when v < 1, W v = ML 

Establishing the upper bound is not difficult, since W v is a limsup set. There 
is no loss of generality in working with the more convenient set W v n [0,1], as 
W v = U kf z Z (W v n [0, 1] + k). Consider the limsup set 

oo oo q 

Win [o,i]= f] U i} B ^>/^T v ~ l l 

N=l q=N p=0 

where B(p/q, e) = {x G [0, 1] : \x — p/q\ < e}, so that dia,m(B(p/q, e)) ^ 2s. Then 

oo q oo 

]T £ diam^b/q, q-^jy < 2 s ^ f-><?+V < oo 

q—1 p—0 q—1 

when s > 2/(v+l). Hence by Lemma, l2~2l a,nd the properties of Hausdorff dimension, 
dimnM^ ^ 2/(v + l) when v > 1. When v < 1, the theorem follows from Dirichlet's 
theorem. 

Establishing the correct lower bound is much harder. Jarnik's lengthy and com- 
plicated proof involved continued fractions and arithmetic arguments. Besicovitch's 
proof was simpler and more geometric and is the basis of regular and ubiquitous 
systems which have turned out to be very effective techniques in determining the 
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Hausdorff dimension of a variety of sets |10L 1131 145j . Indeed ubiquity can imply 
Khintchine's theorem |14| . so ideas developed for the study of the the Hausdorff 
dimension of the null sets also contribute substantially to our understanding of 
this theorem. The Jarmk-Besicovitch theorem has been extended considerably, to 
higher dimensions, hyperbolic space and local fields. Jarmk himself proved the 
Hausdorff measure analogue of Khintchine's theorem for simultaneous Diophantine 
approximation [69j and deduced that the set of points in E™ satisfying 



P 

a 



< 



for infinitely many p € Z", 5 6 N has Hausdorff dimension (n + l)/(v + 1) when 
v ^ 1/n and n otherwise. He also showed that the Hausdorff s-measure of W v at 
the critical dimension is infinite. In view of the inclusions it follows that the 

same holds for E v , R v . The points in W v form an uncountable totally disconnected 
subset of the line and so are in Mandelbrot's picturesque language 'fractal dust', as 
is 03. 

We have seen that the notion of a very well approximable point extends natu- 
rally to systems of real and p-adic linear forms and to hyperbolic space (where they 
can be interpreted in terms of geodesic excursions |33j). The Jarmk-Besicovitch 
theorem and the Hausdorff measure analogue of Khintchine's theorem have been 
established in the real case |41LI421 . the p-adic case [2ll40j . the formal power series 
case |84J and the hyperbolic case 63~l I91j . Other generalisations are to restricted 
sequences |110| . inhomogeneous Diophantine approximation and to small lin- 
ear forms |38| . A further generalisation to 'shrinking targets' has revealed some 
unexpected connections with complex dynamics and ergodic theory [62j. 

We will be interested in the case of a single real linear form in The set of 
f-approximable points a e R™ (see l|3.15|) l will be denoted L v , i.e., 

(4.1) L v — {a: |q • a — p\ < |q|^ for infinitely many qeZ",peZ} 
and by the dual or linear form version of the Jarmk-Besicovitch theorem [22\, 

( , "+ 1 

, . , ,. In— H when v > n, 

(4.2) dim H Lt, = < v+1 

I 11 when v ^ n. 

The n — 1 term arises from the dimension of the resonant hyperplanes 

{x E IT : q • x = p} 
sets while the other represents 'fractal dust' normal to the hyperplanes. 

4.3. Approximation by ratios of Gaussian integers. Recall that the 
Gaussian integers are defined as the set = {pi + ip2 <E C : p\,pi G Z}. These 
form a ring (in fact a unique factorisation domain). The Gaussian rationals are 
defined as the set Q(i) = {a/b + ic/d: a/b,c/d £ Q}. Approximation by Gaussian 
rationals decouples into the independent approximation of the real and imaginary 
part respectively. Here we study the more interesting problem of approximation by 
ratios of Gaussian integers. 

Approximation of complex numbers by ratios of Gaussian integers was studied 
by Hermite and Hurwitz in the 19th century [80, IV, §4] but, unlike in the real case, 
a continued fraction approach did not give the best possible analogue of Dirichlet's 
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theorem. This was obtained in 1925 by Ford |54j, who used additional geomet- 
rical ideas based on the Picard group SL(2, Z[i]). The Gaussian analogues of the 
theorems of Dirichlet, Khintchine, Jarnik and Besicovitch treated below should be 
compared with those corresponding to simultaneous Diophantine approximation in 
the real plane K 2 . 

In 1967 A. Schmidt introduced a theory of regular and dually regular chains 
for continued fractions, to treat approximation problems in complex numbers |111L 
1114] . Schmidt was concerned with the study of complex quadratic irrationals and 
the complex version of Pell's equation and with extensions to groups |112L I113j . 
Our interest, however, is with complex or Gaussian rational analogues of Dirichlet's 
theorem and with the observation that the extended complex plane is the limit set of 
the Picard group. In this connection, Patterson established analogues of Dirichlet's 
theorem for the less general Fuchsian groups |99L 1134) . Later Stratmann and 
Velani obtained versions for Kleinian groups |1_28^ 11291 1135) and so for the Picard 
group. These results can be translated into complex versions of Dirichlet's theorem 
with an undetermined constant. Nevertheless, for completeness, a short geometry 
of numbers proof of the complex version of Dirichlet's theorem is given below. 
Although the constant here is not best possible, the result is all we need. Proofs 
of the complex analogues of Jarmk's theorem on badly approximable numbers and 
the Jarnik-Besicovitch theorem that do not use the hyperbolic space framework 
will be given in ^Jand ^1 below. 

Complex Diophantine approximation has also been investigated from the point 
of view of the distribution of the values of polynomials with real integer coefficients 
but with complex variable z [!§]; for another complex analogue see |5Hj- For the 
rest of this section, p — p\ + ip 2) q = qi + iq2 will denote Gaussian integers with 



Theorem 4.5. Given any z = x + iy e C and N e N, there exist Gaussian 
integers p — p\ + ip2, q — qi + iq2 with < \q\ ^ N such that 



(4.3) 



P 



< 



q\N 



Moreover for infinitely many p,q € Z[i], 



(4.4) 



P 

z 

q 



< 



Proof. The inequality 1)4.311 holds if and only if the inequality 



(4.5) 

holds, i.e., if and only if 



x + iy - 



Pi + iP2 



qi + iq2 



< 



qi + iq2\N 



\{q x x - q 2 y - pi) + i(q 2 x + q x y - p 2 )\ < — 



holds, which is the case if 
(4.6) 



v/2 



max-tlgiX - q 2 y - p t \, \q 2 x + q x y - p 2 \} < -jr- 
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By Minkowski's linear forms theorem, the system of inequalities 



\qix 
\q 2 x- 



Q2V 

Qiy 



,1/2 



AT 



-Pi\ < 

-P2\ < 2 1 /2jV" 1 

\qi\ «S 2~ 1/2 N 
M < 2-V2JV 

has a non-zero solution in integers Pi,P2> <?i> <?2- Hence 114. fifl has a solution with 
\l\ = \li + i<?2| ^ ^V, as claimed. □ 



Since the Gaussian rationals p/q are not required to be on lowest terms, 14.411 
holds infinitely often. For if p/q fails to satisfy H4.4I) . then a fortiori, {np)/{nq), 
where n is a non-zero Gaussian integer, will also fail. And if p/q satisfies 114. 4J1 . only 
a finite number of the fractions (np)/(nq) can also satisfy it. 

This result should be compared to IK-S,7|1 . As in the real case, complex numbers 
for which Dirichlet's theorem cannot be significantly improved are called badly 
approximable and numbers for which it can are called very well approximable. 
More precisely, a complex number z is badly approximable if there exists a constant 
K = K(z) such that for all p,q <E Z[i], q ^ 0, 



P 

z 

<7 



K 



Ford [§4J showed that the complex quadratic irrationals (1 ± iy/S)/2 are the worst 
approximable numbers with K = l/y/3 and thus correspond to the golden ratio 
(y/5 + l)/2 in the real case. Badly approximable numbers have been studied by 
A. Schmidt from the viewpoint of the Markoff spectrum |115| . 
As in the real case, a complex number z is v- approximable if 



(4.7) 



P 

z 

q 



< 



for infinitely many p,q 6 Z(i) and will be called \& -approximable if 



(4- 



<*(M), 



where : [1, oo) — > M + . The set of ^-approximable z will be written W / *( V E'). Thus 
is the set of points z£C for which the inequality l|4.8H holds for infinitely 
many Gaussian rationals p/q. The set of w-approximable complex numbers will be 
written W* . 



4.4. Khintchine's theorem for complex numbers. The complex analogue 
of Khintchine's theorem was proved in 1952 by LeVeque t 85j who combined Khint- 
chine's continued fraction approach with ideas from hyperbolic geometry. In 1976, 
Patterson proved slightly less sharp versions of Khintchine's theorem for Fuchsian 
groups acting on hyperbolic space |99L |l34j. A little later, Sullivan |130j used 
Bianchi groups and some powerful hyperbolic geometry arguments to prove more 
general Khintchine theorems for real and for complex numbers. In the latter case, 
the result includes approximation of complex numbers by ratios a/6 of integers a, b 
from the imaginary quadratic field R(iVd), where d is a squarefree natural number. 
The case d = 1 corresponds to the Picard group and approximation by ratios of 
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Gaussian integers. Stratmann and Velani extended Patterson's results with simi- 
lar minor technical restrictions to Kleinian groups |129L IT35] . These include the 
Bianchi groups and give a less precise and differently formulated version of Sulli- 
van's result. We now state LeVeque's result in our notation. 

Theorem 4.6. Suppose k 2 ^^k) is decreasing. Then the Lebesgue measure of 
W*($f) is null or full according as the sum 

oo 

(4.9) ^2k 3 ^(k) 2 

k=l 

converges or diverges. 

As well as being more general, Sullivan's result is more precise as the growth 
condition for the function ^ is weakened to a 'comparability condition'. Instead 
of the sum, Sullivan and LeVeque use the equivalent integral J x 3 ^(x) 2 dx (their 
integrands are different owing to different forms of H4.8J1 ). It is readily verified that 
W*{^) is invariant under translations by Gaussian integers p — p\ + ip 2 , so that 

(4.10) W*(¥)= (J V(9)+p = (J V*(*)+pi+*p2, 

p€Z[i] pi,P2£Z 

where V*(#) = W*($>) H I 2 and I 2 = [0, l) 2 = {x + iy : < x, y < 1}. We will 
work in the more convenient unit square I 2 and consider the set V*(&). 
Before continuing, we need some more definitions and notation. 

4.5. Resonant sets and balls in I 2 . Let q — q\ + iq 2 € Z[i] \ {0}. The set 

R(q) = R{qi,q2) C I 2 where 

R{q) = l^:p€Z[i\\nI 2 =R{q 1 ,q 2 ) 

(fpm+piqi P2qi-pm\ _ „\ _ t2 

is called a resonant set. This set is the analogue in the complex plane of the set 
{Pi/li '■ ^ Pi < <7i} m the real line. The points in the resonant set form a lattice 
inclined at an angle tan -1 (qi/g 2 ) to the real axis and in which the side length of 
the fundamental region is |g| = (qf + q 2 )^ 1 ^ 2 . Area and congruence considerations 
give that the number of points of R(qi, 92) in I 2 is 

(4.H) #R(q) = #R(qi, 92) = M 2 =q\ + q\- 

The disc 

(4.12) D(p/q; e) = {zeC: \z~p/q\<e}, 
with radius e and centred at p/q where 

P _ Pi + IP2 _ (Pi + jggKgl - ^g2) _ Plgl + P2g2 ^ P2gl -Plg2 

9 9i + «?2 <Zi + g| 9i + ?2 9i + <?! 

has area 7re 2 . The set 

(4.13) B(q,e) = \zel 2 

1 

pel 

can be regarded as a neighbourhood of a resonant set and its measure 
(4-14) \B(q,e)\=ne 2 (\q\ 2 + 0(\q\)). 



P 

z 

q 



< e for some p E Z[i] > = I) D(p/q,e) 
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Now V*(Hf) can be expressed in the form 

oo oo 

(4.15) V*(*)=f) |J |J B(q,*(\q\)) =limsup B(q,V(\q\)). 

N=lk=N\q\=k kl^oo 

Thus V*fif) is a limsup set and for each N = 1,2, . . . has a natural cover 

(4.16) C*r(V*(*)) = {B(q, (|g|) : |g| > TV}. 
Hence for each N = 1,2, ... , the measure of V*(^) satisfies 

oo oo 

rai^E E \B(qMk\)\< E fc2vI/ ( fc ) 2 E L 

fe=A fc^|g|<fe+l fe=JV kt^\q\<k+l 

Now, Xa<M<fe 1 i s the number of lattice points in the closed disc -D(0, k) of radius 
k. By |28j, given any e > 0, 

J2 1 = TTfc 2 + O (fc 12 / 37+£ ) . 

Hence 

(4.17) ^ 1 = 7T(fc + l) 2 + O (fc 1 / 3 ) - TTfc 2 + O (V/ 3 ) = 27Tfc + O (fc 1 / 3 ) . 
fc^|g|<fc+l 

Therefore 

oo 

\v*(*)\ < k 3 ^(k) 2 . 

k=N 

Since iV is arbitrary, the convergence of the series J^k k 3 ^{k) 2 implies that 

\V*(^)\ = |W*(*)| = 0, 

which is the convergence part of the complex analogue of Khintchine's theorem. The 
much more difficult case of divergence requires deeper arguments and the reader is 
referred to Sullivan's bold and highly geometrical paper 130J. 
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5. Badly approximable complex numbers and Jarmk's theorem 

The results of 20J and of [53 on badly approximable points arising with the 
action of Kleinian groups and on bounded geodesies on Riemann surfaces respec- 
tively, specialised to the case of the Picard group, could be translated to give 
Jarmk's theorem for complex numbers badly approximable by ratios of Gaussian 
integers. However, we will give a self-contained proof using an extension to the 
complex numbers of the (a,/3)-game introduced by W. M. Schmidt |117| . In the 
setting of the complex plane, S C C, a,/? € (0, 1) and discs replace intervals. Thus, 
the game begins with player B choosing a disc Bi — {z £ C : | z — b\\ ^ pi}. Next 
A chooses a disc A\ C B\ with radius apt. Then B chooses a disc B2 C A2 with 
radius (3api and so on ad infinitum, such that B n+1 has radius (a(3) n pi for any 
n > 0. 

The discs B\, Ai, B2, A2, ■ ■ . form a nested decreasing sequence of closed sets, 
so there is a unique intersection point, ui say. Player A wins if this point is an 
element of the set S, i.e., if f)j Aj = {u} C S. Otherwise, B wins. A set S is said 
to be (a, /3)-winning if A can win the game for the parameters a and /3 no matter 
how well B plays. If for some a £ (0, 1), A can win the game for any (3 € (0, 1), the 
set S is said to be a-winning. 

Player A benefits from a being small. Indeed, when a gets smaller, A can limit 
the amount of choice B has in the next move. As in the real case, it may be shown 
that if a' < a and S is a-winning, then S is a'-winning. Hence, given S, there is a 
largest value of a for which S is a-winning so we may define 

a*(S) := sup {a € (0, 1) : S is a-winning} . 

In fact, for any S 7^ C, we easily see that a* (S) ^ 5. To see this, note that for 
the parameters a > \ and (5 & (0, 2a — 1), B may ensure that the centres of all the 
Bi are the same. Hence, by choosing the first disc B\ with centre b\ £ S, we have 
the result. 

When S = 25, A wins if she can force u to be badly approximable, i.e., if w is 
in the set 



35 = <^ z e C : 3K > Vp, q S Z[i] 



K 

> 



2 
1\ 



Since 25 is null, this seems unlikely to be the case, but in fact A can win the game 
whenever 2a < 1 + a/3. It immediately follows that 

(5.1) a*(25) = 1/2, 

so in fact A may win this game almost as easily as she could win the game where 
S is the entire complex plane with one point removed. We take some time to prove 

(S3}. 

We have already seen that a*(25) < 1/2, so we will only consider a < 1/2, as 
this will simplify the proof. Let a € (0,1/2] and j3 S (0,1) be fixed. Note that 
7 := 1 + a/3 —2a > 0. We may assume without loss of generality that the radius 
of the initial ball B\ satisfies p\ ^ a/3^. Indeed, otherwise we would let the game 
continue in an arbitrary fashion until reaching a Bj for which pj < a(3^ and then 
take this to be our starting point. The constant K in the definition of 25 will be 
S = j min(p, a 2 p 2 ^), where p = p\. 
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Let t G N be such that a/3j < 2(a/3)' < 7. Let R = (a/3) 
suffices to prove for any n G N that if 



-t/2 



Clearly, it 



(5.2a) 
(5.2b) 
(5.2c) 

then 



> 



gcd(p,q) = 1, 
z G B nt +i, 
0<\q\< R n , 

-. This may be done using induction. 



z - 

1 I 

For n = 0, there is nothing to prove, as l|5.2cjl leaves us no q to consider. 
Hence, we may assume that we have B±, . . . , B(fe_i)t+i such that the above holds 
for O^n^fc— 1. In subsequent play, A thus only needs to worry about | with 
R 1 *- 1 |(jr| < i? fe , as the remaining problematic fractions have been sorted out in 
the preceding steps of the game. 

In fact, there can be at most one such 
P,P', Q, q' G and z, z' G S(fc-i)t+i with 

(5 



Indeed, suppose that there exists 



and 



\q\ 



p_ 

q' 



with |g|,| 9 '| G [R k '\R k ). Then 



+ — o + 2p(a/3) (fe_1)i < 2(5i? 



2-2fc 



k/ 



2 j oi? : 



2-2/c 



7? 



-2A- 



7? 



-2A: 



On the other hand, since gcd(p, q) = gcd(p',q') — 1 and since is a unique 
factorisation domain, 



P_ 

q' 



pq 



p'q 



qq' 



1 



\q\\q'\ 



> R 



-2k 



whenever p ^ p' and q ^ q' . Hence, there can be at most one problematic point. 

Note that we have used the property of the ring of Gaussian integers being a 
unique factorisation domain. For other rings, this may not be the case and stronger 
tools are needed for proving the analogous result. However, for clarity of this expo- 
sition, we will take the simple route and use the unique factorisation property. Note 
also that for the corresponding result in simultaneous Diophantine approximation, 
the module 1? takes the place of the Gaussian integers. Even though the underly- 
ing sets are the same, the difference in algebraic structure prevents the method of 
this proof from working in the case of simultaneous Diophantine approximation. 

As there can be at most one point in B( k _i- )t+1 suitably close to a Gaussian 
rational p/q, we may devise strategies that avoids a disc around this point p/q of a 



suitable radius. It is clear that we need to avoid C 
possibilities in turn. 



6 

w 



). We examine two 



First, consider the case when 



'(fe-i)t+i 



'(fc-i^+i 



> 5R 



2-2k 



> SR 2 - 2k . 
5 



In this case, 



\q\ 
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Figure 3. A avoids C in one move 



Hence, b^-i)t+i t so we are in the situation of figure El As a ^ 1/2, we see 
that A can choose her next disc in such a way that it does not intersect with C. 
This leaves the final case when 



^ SR 



2-2k 



This corresponds to the case when b(f.-i)t+i € C (figureQJ . In this case, we clearly 




B<k-1)H 



Figure 4. A avoids C in t steps 

need to work harder to get an answer, as the right strategy of A is not immediately 
obvious. However, it turns out that picking a fixed direction and moving as far as 
possible in this direction at each turn will cause the disc chosen by B after t steps 
to have empty intersection with C. The following lemma formalises this. 

Lemma 5.1. Let a,/3 € (0, 1) with 7 = 1 + a[3 - 2a > 0. Let t e N be such 
that (a/3)' < 7/2. Suppose that the disc Bk = {bk,Pk) occurs at some stage in an 
(a, 0) -game. Then A can play in such a way that 

B k+t Q{zeC:\z-b k \ >Pfei). 
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Proof. We define a strategy for A in the following way: Suppose that the last 
disc chosen by B was B(c,p) for some p > 0. Choose some c G C with |c| = p(l — a). 
We define a legal move for A by the map 

B(c, p) h-> B(c + c, ap). 

We need to convince ourselves that this move is in fact a legal move. First, we 
note that the radius is the right one. Hence, we need only prove the inclusion 
B(c + c, ap) C B(c, p). But this follows since for any z € B(c + c, ap), 

\z — c\ = \z - c + c — c\ ^ \z — (c + c)| + |c| < ap + p(l — a) = p. 

The strategy of player A will be to use the above move, no matter how B plays the 
game. 

We denote the discs chosen by B by Bk = (bk,Pk), and the discs chosen by A 
by Ak = (a,k,apk). Note that 

(5.3) \a k -b k \ = \b k + c-b k \= p k {l-a). 

Also, from figure we see that 




Figure 5. The point bk+i must be chosen in the shaded area. 

(5.4) \bk+i - dfe| ap k - afipk- 

Hence, by l|5.3^ and I5.4fl . 

\H - b k +i\ = \(bk - a k ) - (bk+i - a k )\ ^ | \b k - a k \ - \b k +i - a k \ \ 

= |/Ofe(l - a) - \b k+1 - a k \ | > \p k (l - a)- ap k + af3p k \ = p k j. 

Continuing as above, at each step choosing c to point in the direction of the first 
one chosen, we obtain, 

\b k+t - h\ > p k ~f. 
But since p k+t = (a/3)V)b < Pfc^, we have for any z e B k+t , 

\z - b k \ ^ | \z - b k+t \ - \b k+t - b k \ \ > p k Z. 
This completes the proof. □ 
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With the induction step done, we have shown that A can play in such a way 
that Ij5.ll) holds for the final intersection point. In fact, we have found an explicit 
lower bound on the constant K for which 115, lj) holds. 

From flQ]) , we can get the Hausdorff dimension of the set 23. By considering 
the amount of choice B has in the game and using this to find a lower bound on the 
s-length of appropriate covers, we obtain a lower bound on the dimension of any 
(a, /3)-winning set. This construction was carried out in considerable generality by 
W. M. Schmidt [I17j, who obtained a lower bound for the Hausdorff dimension 
of (a, /3)-winning sets in real Hilbert spaces. We sketch an approach to obtaining 
the dimension from the number a* (03). Because of the geometric nature of the 
problem, 23 will be regarded as a subset in R 2 . We will compute the usual planar 
Hausdorff dimension here and subsequently. 

To obtain the Hausdorff dimension of 23, we consider the game from B's point 
of view. Assume first without loss of generality that Bg, the first disc chosen, has 
radius 1. At any point of the game, B can choose to direct the game into a number 
of disjoint discs. While there may be a variety of ways in which these discs can be 
chosen, the maximum number N((3) is roughly equal to 1//3 2 , i.e., N(/3) x 1//3 2 . 

We limit B's choice to these N(f3) possible moves and assume that A is playing 
to win the game. This gives us a parametrisation of the sequence of discs chosen by 
B, so that B k = B k {ji, j k ) with the j t € {0, ... , N{j3) - 1} for i = 1, . . . , k. For 
later use, note that the radius of each these discs is p k = (a(3) k . By simultaneously 
considering all the different ways, B may play the game, we obtain a function 

/ : {0, . . . , N((3) - if -> 03, (A fe ) feeN ^ f| B k (X u ...,X k ) = {x(X)}. 

We define the set 23* C 23 to be the range of /. As every number in the interval 
[0, 1] has at least one expansion in base N(f3), we may map this set onto the unit 
interval by the map 

g : 23* -»• [0,1], x(X) h-> O.A1A2 . • . . 

Note that these functions could well be multivalued, but this is of no concern to us. 
All we need is a cover of [0, 1]. We extend this function to subsets of the complex 
plane by defining g(Z) = g(ZC\ < S*) for any Z C C, where by convention £?(0) = 0. 

Now, take some cover C = {C;}/ £ n of 23 with discs of radius p(Ci) = pi. We 
wish to find a lower bound on the s-length of this cover for appropriate s as pi 
becomes smaller. This is where the function g comes in handy. As C covers 23*, 
we see that g(C) covers [0,1]. We let ~p denote the outer Lebesgue measure. By 
sub-additivity, 

00 / 00 \ 

(5.5) E^ C '))>P U»( C ')pi- 

i=i \i=i ) 

Now let lj > be sufficiently small so that any disc of radius Lu(af3) k intersects 
at most two of the discs B k (ji, . . . ,j k ). By |117l Lemma 20], lu = 2/v3 — 1 < 1 
has this property. We define integers 

[ log (2a;- V,) " 
' log(aj9) ■ 



30 



M. MAURICE DODSON AND SIMON KRISTENSEN 



For pi sufficiently small, we see that fc; > and pi < oj(a(3) kl . Hence, the 
disc Ci intersects at most two of the discs B^iji, ■ ■ . , jk t )- As g(Bk l (ji, ■ ■ . , jk t )) is 
clearly an interval of length N((3)- k ' , we have pZ{g{Ci)) ^ 2N(f3) k ' , so by 

oo oo log(Af(/3)) oo log(JV(/3)) 

i ^ ^^(co) < J2 2N W l < 2 ( 2w_1 ) |los(Q/3)l Y,pi los(am ■ 
i=i i=i i=i 

Thus, for s — log(7V(/3))/ |log(a/3)|, the s-length of the cover C is strictly positive, 
so the Hausdorff dimension of 2$* must be greater than this number. 
Now, we fix a e (0, 1/2) and apply the above, 

a- rott lQ S^ 2 2|log/3| 
dlmH(S) > jtoi^ " |loga| + |log/J| 2 

as /3 — > 0. We have thus proved the analogue of Jarnik's theorem for the complex 

numbers: 

Theorem 5.2. The set 93 is thick, i.e., for any disc B C C, 

dim H (Bn<B) = 2. 

6. The complex Jarrrik Besicovitch theorem 

Let W* be the set of w-approximable complex numbers which satisfy the in- 
equality H4.7H for infinitely many Gaussian rationals p/q (recall that p,q £ Z[z]). 
The general Jarmk-Besicovitch theorems in |63L I91| can be specialised to the Pi- 
card group to yield the Hausdorff dimension of W* (see |91l Corollary 2]) but we 
give a more self-contained and direct proof. For convenience we consider 

< — — — for infinitely many — 

Theorem 6.1. 

( 4 

when v ^ 1, 





p 

z 







dimnV v = dining 




when v < 1. 



When v ^ 1, V* is full by the complex form of Khintchine's theorem so the 
second equality holds by (iii) in § 12.21 As usual, the proof for v > 1 falls into two 
parts. First, to obtain the upper bound for dining*, consider the cover c <o{V*) 
given by jlEj with N =1 and &(k) = k~ v ~ l . By (PTTTfl . this has s-length 



DO 



tmv:)) « E E \i\ 2 \\ir~ l ) «E fc2 ^ +1)s E 1 

k=l fc<|g|<k+l k=l fe<|g|<fc+l 



•• E /,;i x 



u3-s(v+l) 
k=l 

for s > 4/(u + 1). It follows that when u > 1, 

(6.1) dim H V; < -4t. 

u + 1 

To obtain the lower bound in the Jarmk-Besicovitch theorem, we use ubiquity 
[HSsl. Let SCI 2 and let p: N -> (0,oo) be a function. Put 

P 
'7 



B(S;q,e) := (zeS: 



< e for some p € Z[i] 
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The set 



K = \jR(q) CC, 



where the union is over non-zero q € Z[i], consists of discrete points and so has 
dimension 0. Let S be any open square in I 2 . By definition, 1Z is ubiquitous in 
S C I 2 = [0, l) 2 with respect to p if 



\jB(S;q,p(N)) 



\S\ 



as TV — > oo. Now by the complex analogue of Dirichlet's theorem, for each JVeN, 

2 



V 

z 



< 



q\N 



9, 



for some p,q € Z[z], 1 < |g| ^ iV 



5'. 



so that 



Consider the set S(N) of z £ S with 'small denominators' 

P 



S(N) = ^zeS: there exist p, q such that 

2 



2 . TV 

< — ,1 < 4 < ; 77 

logiv 



u 



B q 



q\N 



l<|g|<JV/ logJV 

The measure \S(N)\ of S(N) satisfies 

\S(N)\ = \B(q,2/\q\N)\< ]T (\q\N)- 

l<|g|<JV/log N l^\q\<N/logN 



2 M 2 



KA^- 2 (iV/log^) 2 < (logJV)- 



l^|g|<A f /logA' 



as JV — > oo. Choose /o(iV) = 2N~ 2 \ogN. Then since |g| > iV/log7V implies 
p(N)>2/\q\N, 



U 



B(q,p(N))D 



U 



B 9 



AT/ logA r <|g|^Af 



A r /logAf<| 9 |^Af 



\q\N 



s 



in measure as N — > oo and so 1Z is ubiquitous with respect to p(N) = 2 log N/N 2 
for any SCI 2 and so for S = I 2 . But since \& is decreasing, by [45 , 



(6.2) 



dining* ^ dim 1Z + codim 1Z lim sup 



logp(iV) 



where dim is the topological dimension, so that dim 1Z = 0, the codimension 
codimft in C (regarded as R 2 ) of 1Z is 2 and V(N) = N-"' 1 . 

The required result follows on combining the two complementary inequali- 
ties <GlU and ifPJl . 
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In fact, since S was an arbitrary open square, 1Z is locally ubiquitous and it 
has been shown that the local ubiquity of TZ also implies the divergence case of 
Khintchine's theorem |14j . 

Note that even though one might expect Diophantine approximation in C with 
respect to Gaussian integers to be similar to simultaneous Diophantine approxima- 
tion in the real plane M. 2 , this is not at all the case. The analogues of Dirichlet's 
theorem, Khintchine's theorem and the Jarmk-Besicovitch theorem are quite dif- 
ferent in the two cases. Indeed the complex Dirichlet's Theorem and the Jarmk- 
Besicovitch theorem are closer to the real, one-dimensional case. Only the analogue 
of Jarmk's theorem on badly approximable numbers remains unchanged and in this 
case there is a substantial difference in the proofs of the two theorems. 

7. Applications 

The connection between the physical phenomenon of resonance and Diophan- 
tine equations can give rise to the notorious problem of small denominators in which 
solutions to a variety of questions contain denominators that can become arbitrarily 
small. When these small denominators are related to very well approximable points, 
it is sometimes possible to impose appropriate Diophantine conditions which over- 
come the problem by excluding the offending denominators without significantly 
affecting the validity of the solution. The techniques developed in the metrical 
theory of Diophantine approximation lend themselves to this and in particular the 
Jarmk-Besicovitch theorem allows the determination of the Hausdorff dimension 
of the associated exceptional sets. Some examples of problems involving small de- 
nominators and the associated exceptional sets are now discussed. We begin with 
a very simple example. 

7.1. Partial differential equations. Diophantine approximation has been 
applied to the wave equation ( |97j and more recently |51|). as well as to the 
Schrodinger equation |83j . For an extensive treatment of Diophantine problems 
related to partial differential equations, the reader is referred to |104| . 

We will illustrate the applications of Diophantine approximation by Gaussian 
rationals by considering the following innocuous first-order linear complex partial 
differential equation, 

(7.1 a — \-p — ~ = g(z,t), 

at oz 

where z G {x + iy G C : x, y ^ 0}, t ^ and a, f3 are non-zero complex numbers. 
That is, we are studying the partial differential equation on the interior of the set 
defined above under the additional assumption that the functions involved as well 
as all their derivatives may be extended to the whole set. Assume that g(z, t) is 
smooth (i.e., C°°) and can be expressed in the form 

g(z,t) = ^ 9a,b,c,dexp((a + ib)z + (c + id)t) , g a ,b,d,c G C. 

We seek smooth solutions to this equation of the same form, namely 

(7.2) f(z,t)= faAc4^p({a + ib)z + (c + id)t), / Q)M)C G C. 

Thus, we are not just looking for solutions but rather trying to solve the partial 
differential equation subject to boundary conditions. 
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As usual, we solve the problem formally by substituting these two expressions 
into H7.1J1 and identifying coefficients on either side of the equality. Isolating the 
coefficients of /, we get 

1 

(7.3) fa,b,c,d = ~a~ ~~~ : ~~9a,b.c,d- 

£(a + tb) + (c + td) 

We need the coefficients to decay fast enough so that both the series 117. 2JI and its 
derivatives are convergent. Since g is already smooth, the coefficients g a ,b,c,d decay 
rapidly and so are not obstructing this convergence. However, the denominator 
of the fraction may become small and cause the the coefficients f a ,b : c,d to become 
large enough to pose a problem. In order to avoid this, we see that it is certainly 
sufficient for the denominator to be bounded from below by some polynomial in 
a, b, c, d, i.e., we require for some K, v > and for all (a, b, c, d) 6 Z 4 \ {0}, 



(7.4) 



— (a + ib) + (c + id) 
a 



>KmBx{\a\,\b\,\c\,\d\Y 



Since we are only concerned with small denominators, we can assume without loss 
of generality that \a + ib\ x |c + zd|, so that after adjusting K we can drop the 
dependence on c, d on the right-hand side of H7.4J1 and require 

P_P > K 
a q \q\ v+1 ' 

where p = a + ib and q = c + id. Thus we require /3/a to be of complex Diophantine 
type (K, v) for some K, v. The complement of this set is E* — C\ V >\E*, where 

P 

q 

But it can be readily verified by an argument similar to that giving the inclu- 
sion H3.13H . that for each e > 0, 



E* = izeC: for any K > 0, 



K 

< t-. — — for some p, q G Zfil 
\q\ v+1 LJ 



W*^_ C E* C W* 

v v v-\-e ^— v ^— tv v 7 

whence by the properties of Hausdorff dimension given in E 12.2I and the Jarnlk- 
Besicovitch theorem for Gaussian rational approximation (Theorem 16. ljl . 

4 

diniH-E* = — — 
v + 1 

for v ^ 1 and so diniH-E* = lim^oo diniH-E^ = 0. Thus the exceptional set 
associated with the inequality |7.4I failing to hold has Hausdorff dimension zero. 

7.2. The rotation number. The rotation number p(f) is a measure of how 
far 'on average' a continuous, orientation preserving homeomorphism / : S 1 — > S 1 
moves a point round the circle. We will not give the fairly lengthy definition which 
is explained in jJL6j I43L ITU I96j but content ourselves with the observation that the 
rotation number of a rotation r a by an angle 2na, where ^ a < 1, given by 

r a (z) = ze 2ma 

is, naturally enough, a. The rotation number is a nice example of how Diophantine 
properties can arise in analysis as it can be shown that p(f) is irrational if and only 
if / has no periodic points (see |71l Chap. 11,12] or 96, Chap. 1] for more details). 
If p(f) is irrational then for zeS 1 , the closure A of the orbit 

Lu(z) = {f n (z):neN} 
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does not depend on z and either A is perfect and nowhere dense or A = S 1 , 
In the latter case / is transitive and is topological^ conjugate to the rotation 
r p (f) by the rotation number p(f) of /, i.e., there exists an orientation preserving 
homeomorphism tp : S 1 — + S 1 such that 

/ = ( P~ 1 ° »>(/) o <p, 

usually written / ~ r P (f) ■ This can be regarded as obtaining a normal form for / 
and is analogous to diagonalising a matrix. 

Denjoy showed that when / is C 2 and p(f) is irrational, then / is topologically 
conjugate to the rotation by p(f). More subtle aspects arise when additional dif- 
ferentiability conditions are imposed on the conjugation. For example, every C°° 
diffeomorphism / of the the circle is C°° conjugate to a rotation if and only if the 
rotation number p(f ) of / is of Diophantine type |136| . In the analytic case, the 
rotation numbers of real analytic diffeomorphisms form a set lying strictly between 
T> and the set of Bruno numbers f 1371 p. 92]. As in the preceding example, the 
Diophantine condition arises from the denominator 1 — e 27rt P k [ n the coefficients for 
a Fourier series solution of a linearised auxiliary equation in an iterative Newton's 
tangent method argument, modified at each step to retain convergence. In order 
to guarantee convergence of the iterative argument and of the Fourier series, the 
inequalities 



\l-e 2mpk \ ^2 



sin 



fpk-j 



> 2\pk-j\ 



V 2 

where pk—j 6 [0, 2ir), must be set against the very rapid decay of the corresponding 
Fourier coefficient ft- (<C k~ N for any N > 0) in the numerator. It suffices that p 
is of Diophantine type (K, v) for some K > 0, v > 1, since then 

K 



\kp-j\ > j^, i.e. 



J 



k v+l 



for all j/k E Q. 

Now we saw in t!3.2.2l that when v > 1, almost all real numbers are of Diophan- 
tine type (K,v) for some positive K, i.e., the set 

v v= U {a S R: \a-p/q\ > Kq-^ v for each p/q G Q} 

K>0 

is of full Lebesgue measure. Thus the complementary set E v 

E v = P| {a e M: \a-p/q\ < Kq- 2 - y for some p/q € Q} 

K>0 

is null for v > 1 (see M,2.2|) . As in the preceding section, its Hausdorff dimension 
can be determined using the inclusions l|3.13ll and the Jarnlk-Besicovitch theorem 
(Theorem 14.411 

2 

diniH-Eu = dimnT'Fi; = - 

i; + 1 

for v 1. If the rotation number of the smooth circle function / does not lie 
in E — {~\ v> iE v = lim^oo E v (E v decreases as v increases), then / is smoothly 
conjugate to a rotation. The Hausdorff dimension of the exceptional set is 

2 

(7.5) dim H £ = lim = 0, 

v^oo V + 1 

i.e., the complement of V has Hausdorff dimension 0. 
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7.3. The structure of Julia and Fatou sets. Let R(z) — P(z)/Q(z) be 

a rational map on the Riemann sphere C m . A famous result due to Sullivan [90, 
11321 1133j (see also |102j l states that the Fatou set Fr = Coo \ Jr of such a map 
has countably many periodic connected components. These connected components 
were further classified according to the type of periodic cycles they are associated 
with. As the Julia set is the complementary set of the Fatou set, this classification 
also deals with the structure of the Julia set on which the dynamics of Q is chaotic. 

Cycles may be classified according to the value of dR/dz on the points of the 
cycle. By the chain rule, this value remains constant, A say. The cycle is attracting 
(resp. repelling) as |A| < 1 (resp. |A| > 1). Repelling cycles are part of the Julia set 
Jr, so no parts of the Fatou set corresponds to this case. For attracting cycles, an 
associated periodic connected component of the Fatou set is in fact the immediate 
basin of attraction of this cycle, i.e. the union of the connected components of the 
Fatou set containing the points of the cycle in question. 

When |A| = 1, Diophantine properties of A determine the behaviour of the 
dynamics of R and hence the structure of the Julia and Fatou sets. In this case, 
A = exp(27ria). When a is rational, the cycle is said to be parabolic and the 
associated component of the Fatou set is again the immediate basin of attraction. 

In the case when a is irrational, one is interested in the situation when the 
domain of the Fatou set corresponding to the periodic point is a collection of Siegel 
discs associated with the cycle. That is, we are looking for sets on which the 
dynamics are topologically conjugate to a rotation of a disc. Such rotation numbers 
were studied in the preceding section, where the exceptional set associated with the 
failure of this condition was shown to have Hausdorff dimension zero. 

7.4. Linearising diffeomorphisms. Suppose the complex analytic diffeo- 
morphism / : C™ — > C" has a fixed point. Without loss of generality, this can 
be taken to be the origin, so that /(0) = 0. If in a neighbourhood of 0, / is analyt- 
ically (biholomorphically) conjugate to its linear part or Jacobian -D/|o = A say, 
the function / is said to be linearisable. The linearising transformation <fi is given 
by the solution to the functional equation 

/ = (j)' 1 o A o (j>, 

known as Schroder's equation when n = 1. Thus linearisation is similar to conju- 
gating a circle map to a rotation, discussed above in 97.21 And problems of small 
denominators arise when the eigenvalues a\ , . . . , a n of A are close to being resonant 
in the sense that they are close to satisfying the equation 

n 
r=l 

for all j = (ji, . . .,j n ) with j r e N U {0}, r = 1, . . . , n and |j|i = £ r \j r \ ^ 2. 
Linearisation is well understood when n — 1 and the diffeomorphism / : C — > C 
with /(0) = can be linearised when |(-D/|o)| = |/'(0)| 7^ 1. The interesting case 
when I/' (0)| = 1 is closely related via lifts to the conjugacy of a circle map to a 
rotation, discussed in 97.21 above, and necessary and sufficient conditions for the 
linearisation of a diffeomorphism /: C — > C are known [123- On the other hand, 
when n ^ 2, the problem of finding which functions can be linearised is very difficult 
but Siegel's normal forms theorem |123L 1124) gives sufficient conditions on Df\ 
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for the existence of a linearising transformation <fi. The point (a%, . . . , a n ) in C™ is 
said to be of multiplicative type (K,v) p. 191] if 



(7.6) 



r=l 



for all j e (NU {0})™ with |j|i ^ 2. Siegel showed that if the vector (ai, . . . , a„) of 
eigenvalues of Df\ is of multiplicative type (K, v) for some K > and u > 0, then 
/ can be linearised in a neighbourhood (further details are in fTl Kill 1931 ). To stop 
Siegel's condition being too restrictive, one chooses v > (n — l)/2, since then the 
set of points of multiplicative type (K, v) has full measure for any K > 0. However, 
the neighbourhood of linearisation decreases as v increases (it also depends on K 
but less significantly) and so we do not want v to be too large. 

Let £ v denote the exceptional set of points in C" (regarded as IR 2 ™) which for 
a given exponent v, fail to be of multiplicative type (K, v) for any K > and so 
fail to satisfy the conditions of Siegel's theorem and suppose if v > (n — l)/2. Then 
£ v is null and its Hausdorff dimension is given by 14.21 ; namely, 

71+1 

(7.7) dim£„ = 2(n- 1) + 



' + 1 

This result is established by means of an exponential map which preserves the 
Hausdorff dimension and allows £ v to be replaced by a set involving a more general 
kind of additive type (see next section) with a simpler structure |46j . 

7.5. Lyapunov stability of vector fields. Consider the differential equation 



(7.8) z = Az + Q(z)eC, 

where A is a n x n complex matrix and the holomorphic functions Q : C" — > C™ 
and dQk/dzj vanish at the origin 0. The obvious solution zo(t) = is said to be 
future (resp. past) stable if points near remain there under evolution by 117. 8J1 
to the future (resp. past). More precisely, the solution z(t) is future (resp. past) 
stable if for every neighbourhood TV of 0, there exists a subneighbourhood N' with 
e N' C N such that z(0) G N' guarantees that z(i) € N for all* > (resp. t < 0). 
By a well known theorem of Lyapunov |93j , the solution is future stable if the real 
parts of the eigenvalues of A are at most and past stable if the real parts are at 
least 0. Thus for the solution z(t) to be future and past stable or simply stable, the 
eigenvalues must have zero real part and so must be purely imaginary. The stability 
of the solution z(t) is determined by a remarkable theorem due to Caratheodory and 
Cartan |93j which asserts that stability is equivalent to A being diagonalisable with 
purely imaginary eigenvalues and the vector field being holomorphically linearisable 
in a neighbourhood of the origin. By Siegel's normal form theorem |123LlT24] . also 
used in the preceding section, this last condition holds if the eigenvalues "fk = iXk, 
Afc G M, satisfy for some K>0, v > n — 1, 



(7.9) 



^fe ~ y ' A r j r 

r=l 



for all j G (NU {0})" with |j|i ^ 2 (note that this is the additive form of fTffl ). 
The complement of this set of points A e W 1 of additive type is the set of a E R" 
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such that for any K > 0, 
(7.10) 



ttfc-^ar jr 

r=l 



<*ijir 



for some j € Z". This set, E v say, is related to and has the same metrical character 
as the set 

(7.11) L v — {a: |q ■ a\ < |q|^ for infinitely many q€P} 
and in fact by an argument similar to that giving H3.13J1 . for any e > 0, 

L v j re C E v C L v 

(see |16t Sect. 7.5.2]). This inclusion implies that the two sets L v and E v have the 
same Hausdorff dimension. The set L v is also related to L v and roughly speaking, 
has one degree of freedom less. The Hausdorff dimension of L v is a special case of 
an 'absolute value' analogue, proved by Dickinson 38j, of the general form of the 
Jarnik-Besicovitch theorem: 

(7.12) dim£„ = dimE v =n-l + 



+ 1 

when v > n— 1 (see also |47| ): note that L v = R™ otherwise. Thus the exceptional 
set of eigenvalues for which the solution to l|7.8|) cannot be shown to be stable has 
Hausdorff dimension 0. 

7.6. Kolmogorov-Arnol'd-Moser theory. The stability of the solar system 
is one of the oldest problems in mechanics 94 . It is of course a special case of the 
N body problem of understanding the motion of N point masses subject only to 
gravitational attraction, with all other forces neglected. When N = 2, the solution 
is well known and the periodic solutions in which the bodies move in an ellipse 
about their centre of mass persist forever. For N ^ 3, however, the situation is 
extraordinarily complicated and is far from being fully understood, even for solar 
systems where the mass win of the sun is much greater than the masses of the 
n = N — 1 planets. If, as a first approximation, the centre of mass of the system 
is assumed to coincide with that of the sun and if the gravitational interactions 
between the planets and other effects are neglected, the system decouples into n 
two-body problems, in which each planet describes an elliptical orbit around the 
sun, with period Tj say and frequency uij — 2ir/Tj, j = 1, . . . , n. 

For each vector lo = (u>i, . . . , w n ) of frequencies in the n-dimensional torus 
T" = S 1 x • • ■ x S 1 , the map ip w : M -> T n given by 

tp w (t) = <p u (0) + 

is a quasi-periodic flow on the torus. 

The case n = 1 corresponds to uniform motion around a circle and so is peri- 
odic. When the frequencies are all rational, the flow is periodic. If the frequencies 
are not all rational, then by Kronecker's theorem [57 j5 the flow winds round the 
torus, densely filling a subspace of dimension given by the number of rationally 
independent frequencies. Thus when the frequencies are independent, the closure 
of is the torus T n and solutions will persist for ever. Gravitational interac- 

tions between the planets are represented by a small perturbation of the original 
Hamiltonian describing the system. Stability then reduces to the solutions of the 
perturbed Hamiltonian system continuing to wind round a perturbed invariant 
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torus. Of course this model is idealised and takes no account of the final fate of the 
universe. 

Details of the history of the solution to this problem are in [93, Chap. 1]. 
Siegel's success in overcoming the related 'small denominator' problem in the lin- 
earisation of complex diffeomorphisms (see H7.4I M7.5I1 was followed by Kolmogorov's 
conjecture that quasi-periodic solutions for a perturbed analytic Hamiltonian sys- 
tem not only existed but were relatively abundant in the sense that they formed a 
complicated Cantor type set of positive Lebesgue measure |81j . This was proved 
completely in 1962 by Arnol'd pj] and independently Moser proved an analogous 
result for sufficiently smooth 'twist' maps [92, 125J. The results imply that for 
planets very much smaller than the sun and for the majority of initial conditions 
in which the orbits are close to co-planar circles, distances between the bodies will 
remain perpetually bounded, i.e., the planets will never collide, escape or fall into 
the sun. Further details can be found in [8 J and 36j. 

The differentiability and Diophantine conditions were relaxed substantially by 
Riissmann in [ToTl ITok Il09j (see also H Sect. 6.3], gE], [Ml Chap. 1], [103J). 
Another approach is to use 'averaging' methods pp; this can involve Diophantine 
approximation on manifolds |44| . see also 109J. 

As in the above examples, it turns out that in order to ensure convergence of 
a Fourier series and an infinite dimensional extension of Newton's iterative tangent 
method, the frequencies lj — (u>i, . . . , u> n ) must satisfy a Diophantine condition, 
which in this case is 

(7.13) |q- u\ = {qtU! + ■■■ + q n u> n \ > K\q\^ v , 

for some positive constants K — K{u) and v — v(ui) for all non-zero q = (qi, . . . , q n ) 
vectors in Z™. The exponent v is subject to two conflicting requirements. It should 
be large enough (v > n — 1) to ensure that the Diophantine condition above is 
not too restrictive, but small enough to ensure that the perturbation has physical 
significance and that the stability is robust. The proof breaks down when the 
frequencies lie in the complementary exceptional set E v , say, of frequencies which 
are close to resonance in the sense that, given any K > 0, there exists a q e Z n 
such that 

|q-a| <K\q\r- 

The set E v is related to the set L v discussed in the preceding section and allows us 
to deduce that 

diniH-E^ = dimn^ti = n — 1 H 

v + 1 

when v > n — 1. 

Thus Hausdorff dimension plays its part in other branches of mathematics and 
mechanics, as well as in number theory. It has even been of use to Mandelbrot in 
arousing the interest of mathematicians in his profoundly original and imaginative 
ideas 29 Chapter 2]. 
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